Modulated superconductivity due to vacancy and magnetic order in A x Fe 2 - x /2^2 [A =Cs, K, 

(Tl,Rb), (T1,K)] iron-selenide superconductors 

Tanmoy Das 1 , and A. V. Balatsky 1,2 
Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA. 
2 Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA. 

(Dated: June 17, 2011) 

We present a calculation of a 'modulated' superconducting state in iron-selenide superconductors. The zero- 
momentum d— wave pairing breaks the translational symmetry of the conventional BaFe2Se2-like crystal of 
IA/mmm space group. This pairing state becomes compatible when the Fe vacancies form an ordered state 
and the crystal symmetry changes to a low-temperature 14/ m one. For the specific case of an incommensurate 
vacancy order at Q v — (i, |) in Ko.82(2)Fei.626(3)Se2, we find that it induces a block checkerboard antifer- 
romagnetic phase at wavevector Q m = 4Q„. The coexistence of vacancy order and magnetic order leads to a 
reconstructed ground state which naturally couples to the d— wave superconductivity in a uniform phase in what 
we propose will be a general coupling for all iron-selenide superconductors. Our results agree with numerous 
experimental data available to date. We thus suggest that the incommensurability leads to uniform coexistence 
of multiple phases as a viable alternative to a nanoscale phase separation in high-r c superconductors and play 
an important role in the enhancement of superconductivity. 



PACS numbers: 74.81.-g,75.10.-b,74.20.Rp,74.70.Xa 

Introduction— The layered iron-selenide compounds 
A x ¥e, 2 _ x / 2 St 2 [A =Cs, K, (Tl,Rb), (T1,K)] have been 
found recently to become superconducting (SC) with highest 
T c ~ 32K when a suitable value of the Fe vacancy concentra- 
tion x ~ 0.87 is achieved.] 1 1 This fascinating tunability of the 
SC property with Fe vacancy has opened up a new platform 
to study the mechanism of unconventional superconductivity. 
Most interestingly, the onset of vacancy order is accompanied 
by several interesting changes of the structural, electronic, 
magnetic and SC properties. In particular, numerous experi- 
ments have found that the superconductivity turns on exactly 
when the randomly placed Fe vacancies also form an ordered 
state in addition to a structural phase transition. 1 2 1 Simulta- 
neously, the value of the magnetic moment is dramatically 
enhanced on each Fe sites. 

We show here that in these materials the vacancy or- 
der induces the magnetic order, and that both phases uni- 
formly coexist with ci-wave superconductivity to lead to a 
'modulated superconductor'. The induction of a modulated 
order due to the incommensurability of another, compet- 
ing, order is a common phenomenon in unconventional su- 
perconductors. Historically the incommensurate SDW or- 
der in cuprates was first found in neutron scattering experi- 
ment La! 6 _ 2; Ndo.4Sr 2 .Cu04 with wavevector Q m = (it ± 
5, 7r)/(7r,7r±(5),[3 | similar to the case ofCr.[4| Subsequently, 
an induced CDW order at Q c = (2d, 0)/(0, 25) is also mea- 
sured in scanning tunneling microscopy (STM).||5] |6) Coex- 
istence of SDW and CDW produces a 'stripe' -like pattern in 
cuprates which is observed by neutron scattering and ARPES 
in La-based cuprate, STM in Bi2201 and quantum oscilla- 
tions, Nernst signal and Seeback measurements in YBCO 
families. 042) Similarly, in iron-pnictides (BaFe2As2) and 
chalcogenides (FeTei-^Se^), the SDW with a modulation 
Q m = (tt, 0)/(0, tt) has been argued to lead to a CDW 
with Q c — 2Qc. lfl0l The phenomenology also extends to 



the 'hidden-order' state of heavy fermion URu 2 Si2 which has 
an incommensurate modulation at Q/j = (0.6-7T, 0)/(0, 0.67r) 
and the formation of the CDW at Q c = 2Qh has been 
predicted- full It is thus timely and important to investigate 
the nature of the density wave order of the conventional SC 
quasiparticles, and their coupling to the vacancy order in Fe- 
selenide as another example of same basic principle of modu- 
lated states high-T c superconductivity. 

Fig. 1(a) gives the experimental phase diagram! 1 J as a 
function of temperature and Fe vacancy for A E Fe 2 _. I ;/2Se2. 
Our calculations consistently demonstrate that the vacancy 
order in Fe-selenide at a superlattice modulation of Q„ = 
(|j | ) El El breaks the translational symmetry of the tetrag- 
onal crystal in such a way that it induces structural, magnetic 
and unconventional SC phase transition exactly at the same 
value of x ~ 0.87. At the vacancy order transition tempera- 
ture T v ~578K in K .g2(2)Fe 1 .626(3)Se2[2|, the crystal struc- 
ture undergoes a transition from the high-temperature tetrag- 
onal BaFe2As2-like structure of IA/mmm space group to a 
lower symmetry tetragonal structure of 14/ m where the lat- 
ter can be understood as a y/b x \/5 x 1 superlattice of the 
former, see Fig. 1(b) Subsequently, a block checkerboard like 
antiferromagnetic (AFM) order stabilized with significantly 
enhanced moments on the Fe sites [about 3.3/is] below a 
Neel temperature Tn sa 559K< T v with an order vector 
Qm = (|i |)-[2| Finally, the superconductivity appears be- 
low T c ~ 32K. Later on this overall phenomena has been 
observed to be a common feature in all iron-selenide chalco- 
genides studied so far.lffl frZl 

Vacancy order and magnetic order— We concentrate here 
at a particular case of K .82(2)F e i.626(3)Se2l21 at which the 
vacancy superstructure and the magnetic structure have prop- 
agation vectors Q„ = (\, §) and Q m = (§,§). The 
property of the magnetic order Q m = 4Q V motivates us to 
study a modulated order model in which the vacancy order 
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FIG. 1. (Color online) (a) Experimental phase diagram of A iE Fe2_x/2Se2 compound taken from Ref. (T] is also shown to be a common phase 
diagram for this class of materials A c Fe 2 _ a; /2Se2 in Ref. 1121 . At high temperature, the system is in the Ba2Fe2As2-like tetragonal structure 
of 14/mmm space group at all Fe vacancies. A tetragonal structure with 14/ m unit cell is observed to stabilize at low temperature in the 
v5 x v5 x 1 superlattice. The metal-insulator transition is observed around x ~ 0.865 in numerous experiments albeit a smooth behavior 
of the Neel temperature Tat (diamond) throughout the phase diagram. The open diamond symbol marks the transition temperature at which 
disorder to order state of the Fe vacancy is measured. Solid squares represent T c ~ 30K which seem to be symmetric around x = 0.8 at 
which full occupancy of Fe2 sublattices and vacant Fel sites is realized, (b) Common crystal and magnetic structure of AsF&iSes in the real 
space tetragonal 14/ m unit cell (top view), only Fe atoms are shown.fT1l2l ll2l Magenta and cyan color circles give spin-up and spin-down 
configurations, while black and black dumbbells represent positive and negative phase of the pairing symmetry. Boxes of different colors and 
sizes stand for different unit cells (text), while yellow box is the true unit cell for this system at low-temperature. 



is the driving instability which induces the magnetic order 
at Tn < T v . Vacancy order will naturally develop a CDW 
at 2Q V . Neutron and x-ray study in a Cs intercalated com- 
pound CSj,Ke 2 _,j./ 2 Se 2 showed the presence of the superstruc- 
ture reflection at a Q = (|, |)|T4'] which is equal to 2Q V 
(antiphase). Finally, the vacancy order clearly breaks the rota- 
tional C\ symmetry of the crystal, leading to a global nematic 
order parameter. The unidirectional nematic order breaks the 
t2 g symmetry of the Fe d orbitals (which contribute mostly to 
the FS) and thus an orbital order around each vacancy can also 
be expected. We assume the orbital order will have modula- 
tion at 3Q V . Such an induction process is reversible as long as 
the modulation vector is incommensurate which we assume to 
be the case. It is interesting to note that such possible nematic 
order is consistent with claims of nematicity in cuprates,|7), 
pnictideEOl, URu 2 Si 2 |[n] Q2I and Sr 3 Ru 2 7 |[l6), although 
the details of the ordering are material sensitive. 

The coupling between the modulated order with d— wave 
superconductivity are strongly constrained by symmetry 
and momentum conservations and belong to a point-group- 
symmetry. Due to C4 symmetry breaking, s-wave and d- 
wave representations always mix, but d x i_ y i and d xy pair- 
ing remain clearly distinguishable. Q El Therefore, the order 
parameter can be defined globally via the electronic momen- 



tum distributions as faj = J2k^( c k' a ^a,/S c k^), where d^ 
is the structure factor for various d— waves and S Q ( g = <5 Q ( g 
or Ca,B f° r scaler and vector order, respectively. Decou- 
pling the global parameter into the local orders, one can ac- 
quire well defined directions as fa = ^ k £r (cj c+ Q. (7 Cka-) for 

scalar or fa = X)k ct( c 1c+q t Ck ^ ^ or s P' n orc ^ ers - Simi- 
larly for superconductivity, we obtain modulated orders ip^ = 
Sk( c ik+Q -t c -k^' while the zero-momentum SC order pa- 
rameter is ipo = X)k( c iLt c -H)- I n tn e present case, we re- 
strict ourselves to case with i = l — 4:i = l for vacancy 
order, i = 2 for charge order, i — 3 for orbital order and 
i = 4 for magnetic order. Keeping only the linear terms in the 
global modulated order parameter <fi n in the tetragonal crystal 
environment, we can write all the coupling terms Q 

4 4 
n 

+ J2 [A^i# + Af(^Vo^ + ^i^o#)... + c.c] .(1) 
»=i 

The coupling constants Ah between vacancy order with oth- 
ers are possible only if the corresponding modulation follows 
Qi = nQi relation (Qi = Q„). We will comment on the 
superconductivity below. 
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In Fig. 2, we show the evolution of the FS reconstruction 
and the nature of the gap opening in the normal state. Ab- 
sorbing the coupling constant in the order parameters, we de- 
fine the mean-field gaps as tyj (= Aj</>,-). The eigenvector is 



*S = (4, CT C l+ Ql ,a C k+Q 2 ,<r> -)' where C k+ Ql , ff is the cre- 
ation operator for the Qith order and so on. We construct the 

two band model Hamiltonian in 2 Fe unit cell in the normal 

state [i.e. when SC order parameters ipo and ip vanish] as 
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Bare dispersion, £fc, matrix arises from two tight-binding 
bands of strongly hybridized t 2g — orbitals present near the 
Fermi-level. [ 18 ] In such incommensurate case, the FS recon- 
struction is determined by multiple translations of the bare 
band structure £ k by ±nQ v , where the largest value of the 
integer n is determined by Q t , = (2ir, 2tt) in 2D which is 
10 in our case. Each component of the translated disper- 
sion is denoted by £k+«Q- As in the case of non-SC Cr,[4| 
cuprate|17| or other SC materials iflOl fTTl [T5l [l6l incommen- 
surability causes hierarchy of gap openings of order A^f ~ 
[V» m lW m - x + V^ m /W m - X ) [W is bare band width, A 
should not be confused with SC gap] at band crossings be- 
tween ££ +nQ and £k+ (n±m)Q .03D As long as Vi « W, the 
FS is well described by including the lowest-order gap only 
and we will neglect all the matrix elements with m > 1 and 
the momentum dependence in the potentials Vi. 

We have used the eigenstates and eigenvectors of Eq. 2 to 
compute the band dispersion along high-symmetry lines and 
the corresponding FS in the normal state, as shown in Fig. 2. 
The results are presented for various Vi, while fixing the num- 
ber of quasiparticles on the FS. We also present the real part 
of the non-interacting static susceptibility in the 2D g-space 
ignoring matrix-element effects. 

To understand the evolution of the modulated order due to 
turning on of individual interactions, it is useful to examine 
them first separately. In the nonordered state (all Vi = 0), 
the FS consists of two electron-like closed orbits centering 
M— point in the 2 Fe unit cell, see Figs. 2(al,a2). The corre- 
sponding susceptibility in Fig. 2(a3) shows nesting peaks or 
streaks at all nQ v s while the intensity dominates at Q v and 
Qm = iQv With turning on of V v = — 0.5eV, the vacancy 
order is illuminated and the FS reconstruction becomes com- 
plicated, see Figs. 3bl,b2,b3. While both electron pockets are 
still present in closed orbits, but their size and the associated 
spectral weight are reduced. The missing spectral weight is 
shared by a newly developed electron-pocket and several open 
orbit FS pieces. The phase of the interaction V v does not have 
significant effect on modulated pattern (not shown). 

We now study the magnetic ordering in Fig. 2(c). To con- 



nect the magnetic order parameter V m with experiments we 
will investigate local moments. Neutron diffraction has de- 
termined a large magnetic moment block checkerboard AFM 
in the range of 2.3-3.3/is/Fe for all known iron-selenide 
superconductors. l2l [T2ll Within our mean-field approach, we 
find that a value of V m = 20meV produces maximal lo- 
cal moments of 2/Lts(S' z ) max « 0.12/is in 2 Fe unit cell 
which is smaller than seen in experiments. The value of com- 
puted magnetic moments can be enhanced by various modi- 
fications of our calculation, e.g. when all d— orbitals are in- 
cluded. With all interactions included, a drastic band folding 
and FS reconstruction are observed in Figs. 2(d). The cor- 
responding susceptibility continues to exhibit dominant nest- 
ing along the vacancy order vector. First principle calcula- 
tion have demonstrated drastic FS reconstruction driven by 
the Fe vacancy and AFM order. 1 20 1 Such a FS reconstruction 
is consistent with observations such as transmission electron 
microscopy (TEM),lfT3l NMR||2TJ, Neutron diffraction (2) and 
optical studies. [22 1 

Vacancy order and d-wave pairing. — As superconductivity 
turns on, the coupling between the competing order to SC Af 
in Eq. 1 are allowed only when Qi = iQ s where i = 1, 2... 
(Q s is the SC ordering vector in the pair-density channel). In- 
terestingly, the zero-momentum condensate Q s = is always 
a complimentary component of the SC state, whereas the finite 
momentum components are subject to the higher order cou- 
pling strength and depend on the nature of the driving insta- 
bility vector Qi. Note that Xf alone generates a Fulde-Ferrell- 
Larkin-Ovchinnikov state. [ 23 1 Absorbing the coupling con- 
stants Af in the order parameter, we define the SC gaps as 
Ai = Xftpi- Similar to Eq. 2, we can construct the SC gap 
matrix H sc = (A k A k+Ql A k+ q 2 ...;c.c). In the Nambu 
space defined by *t(k) = (*J(k); ^(-k)), we get the 
modulated SC Hamiltonian as 



H = * T 



H (k) tf sc (k) 
H* sc (-k) ff (-k) 



(3) 



In the rest of the paper we discussion how the vacancy 
order play an important role to the onset of d— wave super- 
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FIG. 2. (Color online) The band dispersion along high-symmetry lines [(al),(bl),(cl),(dl)], corresponding FSs [(a2),(b2),(c2),(d2)] and bare 
susceptibility [(a3),(b3),(c3),(d3)] maps at zero energy are presented for four combinations of interactions: (a)— ¥ all Vi = 0, (b)— > only 
V v = -0.5eV, (c)-> only V m = 0.02eV and (d)-> V v = -0.5, V 2 = 0.2, V" 3 = -0.1, V m = O.OleV. FS and susceptibility are plotted only 
in the first quadrant of the BZ which are symmetric in the other quadrants. Note that opposite phases between consecutive interactions give 
maximal weight on main bands. Note that the electronic structure consists of 20 split bands while the spectral weight is largest in the main 
band if coherence factors of the modulated order are included. Here we have not included the structure factor of the modulated order which 
may further enhance the distribution of the spectral weight among the modulated bands. 



conductivity. It has been argued that the d— wave order in 
these materials without a hole-pocket at V in the conven- 
tional BaFe2As2-like crystal structure breaks the translational 
symmetry. 1 18, 24 1 To elucidate this, we recall first the exam- 
ple of BaFe2As2 pnictide superconductor which hosts both 
hole and electron pockets at T and M points, respectively, in 
the 2 Fe unit cell and how s ± -pairing is compatible with the 
corresponding crystal symmetry. We illustrate this with the 
help of the unitary transformation from 1 Fe unit cell to the 



2 Fe unit cell where one expects, by construction, to retain 
all observables such as nodeless SC gap(s) and FS nesting to 
be compatible with the actual crystal symmetry irrespective of 
the conventional one used for calculations. 

We define a quantity in the 2 Fe unit cell by tilde. The 
unitary transformation consists of 45° rotation of the crystal 
with lattice constant a —> \/2a which gives k x / y = (k x ± 
fcy)/2.[25| In doing so, the s ± -pairing in the 2 Fe unit cell 
becomes a s^-like symmetry in the 1 Fe unit cell in which 
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FIG. 3. (Color online) Schematic behavior of the FS and pairing symmetry in 1 Fe per unit cell and 2 Fe per unit cell Brillouin zone for 
BaFe2As2 are shown in (a) and (b), respectively. The corresponding real-space view of the gap symmetry is depicted in (c) in which the SC 
phase is isotropic and nodeless on each Fe atoms but changes phase between Fel and Fe2 sublattices. Black to white color scale depicts the 
pairing symmetry from negative to positive in momentum space while black and white dumbbells represent positive and negative phase of the 
pairing symmetry in real space (exponential form of cosine and sine functions). Similar results for AFe2Se2 compounds are given in (d)-(f). 
Although in momentum space FSs are isotropic and nodeless, but in real space the gap nodes pass through both of the Fe atoms. The phases 
of the pairing symmetries are s xy — 2cos (k x a) cos (k y a), s ± /d x 2_ y 2 = cos (k x a) ± cos (k y a), and d' xy — 2sin (k x a) sin (k y a). 



two inequivalent Fe sites posses opposite phases in real space, 
see Fig. 3(a-c). In momentum space, the SC gap symmetry is 
compatible with the FS topology. In real space, the gap phase 
obeys proper translational symmetry in the true 2 Fe unit cell, 
but not in the 1 Fe unit cell (conventional). This seems to pose 
a problem that is in fact an artificial issue. Proper unit cell 
contains 2 inequivalent Fe sites; 1 Fe unit cell is an artificial 
one constructed for convenience of calculations assuming 2 
inequivalent Fe to be the same. 

We now look at the SC gap transformation. As obvious in 
pnictide, the sign change of the SC gap in s ± -pairing occurs 
between Fel and Fe2 atoms, see Fig. 3(c). Therefore, when 
the center Fe2 atom is removed from the lattice (when vacan- 
cies are present and ordered) in A x ¥s 2 -x/2^2, one can ex- 
pect a different SC pairing symmetry to obtain sign change 



in the bulk SC order parameter. All of the Fel atoms are 
now equivalent, the sign change of gap in the pairing func- 
tion should be carried by all atoms, i.e. a d— wave or p— wave 
state can arise. lnA x Fe 2 -x/2S s 2 compounds when T c is high- 
est at x ~ 0.8, a d— wave gap is a most natural pairing state 
to arise. [18 1 NMR studies also indicate the presence of spin- 
singlet pairing channel in the K intercalated samples with dif- 
ferent composition. |2T) The d— wave pairing which gives rise 
to nodeless and nearly isotropic SC gap on the FSs, is also 
consistent with numerous observables. lfT8l 

On the other hand the d— wave SC gap in A; C Fe2_ :r /2Se2 
compounds has an additional complication. While a d x 2_ y 2- 
wave gap is compatible in both momentum and real space of 
the conventional 1 Fe unit cell, it is inconsistent when the pair- 
ing symmetry is transformed to the 2 Fe unit cell. For exam- 
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pie, in Fig. 3(e), when a k — (— 7T, — 7r) point is shifted by 
the reciprocal vector of (2ir/a, 0) (assuming a square lattice), 
the FS topology is restored but the SC order maps on the site 
with an opposite sign of the gap function. We propose that 
this problem can be cured by doubling the Brillouin zone in 
the momentum space. One way of doing it is by removing the 
Fe2 atom from the center of unit cell in Fig. 3(f). Therefore, a 
d— wave gap is possible in these compounds when the normal 
state vacancy order helps restore the translational symmetry 
of the crystal. 

Iron-selenide materials host a SC dome-like feature center- 
ing at x « 0.8 as long as vacancy order persists, see Fig. 1(a). 
It has been reported that at a; = 0.8 all Fe2 atoms from the 
sublattice are completely removed, [2]. This observation is 
consistent with the symmetry expected for the d— wave pair- 
ing. Such a pattern of competing order is forbidden in the 2 
Fe 122-pnictide IA/mmm structurel 12] and the system stabi- 
lizes to a low-temperature tetragonal 14 /m one, as indicated 
by the yellow box in Fig. 1(b). 

In conclusion, we present a microscopic model to show that 
the onset of vacancy order leads to many significant modifi- 
cations of electronic, magnetic and superconducting proper- 
ties of iron-selenide superconductors. A block checkerboard 
AFM state with moments as large as 3.3/is is realized ex- 
perimentally and is believed to be induced by the vacancy or- 
der at Q m = 4Q„. Such a modulated ground state compen- 
sates for the broken translational symmetry of the d— wave 
superconductivity! 18 1 and therefore allows for the state with 
combined superconducting and structural modulations to de- 
velop. In this state superconducting gap is nodeless and nearly 
isotropic on the FSs, as observed in numerous experiments. 

The results follow the same route as in other high-X^ super- 
conductors where a modulated ground state of multiple com- 
peting phases promotes and coexists with superconductivity in 
a uniform phase rather than undergoing any nanoscale phase 
separation. The modulated superconductivity has also been 
seen in cuprates.ll7l [T71 pnictides.lflOl URu 2 Si 2 |[ni S3 and 
Sr^Ru9Q7 lfT6l . although the nature of competing phases may 
vary across these materials. 
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